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ABSTRACT 



Context. The region G24.78+0.08, which is associated with a cluster of high-mass young stellar objects in different evolutionary 
stages, is one of the best laboratories to investigate massive star-formation. 

Aims. We aim to image the molecular outflows towards G24. 78-1-0.08 at high-angular resolution using SiO emission, which is consid- 
ered the classical tracer of protostellar jets. In this way we study the mass loss process in which we previously detected a hypercompact 
ionised region, as well as rotation and infall signatures. 

Methods. We performed SiO observations with the VLA interferometer in the / = 1-0 v=0 transition and with the SMA array in the 
5-4 transition. A complementary IRAM 30-m single-dish survey in the (2-1), (3-2), (5-4), and (6-5) SiO lines was also carried out. 
Results. Two collimated SiO high- velocity (up to 25 km s"' w.r.t. the systemic velocity) outflows driven by the A2 and C millimeter 
continuum massive cores have been imaged. On the other hand, we detected no SiO outflow driven by the young stellar objects in 
more evolved evolutionary phases that are associated with ultracompact (B) or hypercompact (AI) Hii regions. The A2 outflow has 
also been traced using H2S. The LVG analysis of the SiO emission reveals high-density gas (lO^-lO'* cm"^), with well constrained 
SiO column densities (0.5-1 10'^ cm"^). The driving source of the A2 outflow is associated with typical hot core tracers such as 
CH3OCHO (methyl formate), C2H3CN (vinyl cyanide), HCC'^'CN (cyanoacetilene), and (CH3)2CO (acetone). 
Conclusions. The driving source of the main SiO outflow in G24 has an estimated luminosity of a few 10* Lq (typical of a late 
0-type star) and is embedded in the 1.3 mm continuum core A2, which in turn is located at the centre of a hot core that rotates on a 
plane perpendicular to the outflow main axis. The present SiO images support a scenario similar to the low-mass case for massive star 
formation, where jets that are clearly traced by SiO emission, create outflows of swept-up ambient gas usually traced by CO. 

Key words. ISM: individual objects: G24.78+0.08 — ISM: molecules — stars: formation 
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1. Introduction 

Two main theoretical scenarios, based on accretion, are proposed 
to explain the formation of O-B type stars: (i) the core accretion 
model (McKee & Tan 2002, 2003), where massive stars form 
from massive cores, and (ii) the competitive accretion model 
(Bonnell et al. 2007), where a molecular cloud fragments into 
low-mass cores, which form stars that compete to accrete mass 
from a common gas reservoir. Both models predict the existence 
of accretion disks around massive young stellar objects (YSOs), 
and the presence of jets driving molecular outflows. The core ac- 
cretion model is a scaled-up scenario of low-mass star formation. 
The competitive accretion model suggests that massive stars al- 
ways form in densely clustered environments and that disks and 
collimated jets are perturbated by interaction with stellar com- 
panions. Observation of YSOs with disk/jet systems, and of their 
properties, would help to distinguish between models. 

The region G24.78-f0.08 (hereafter G24), located at 7.7 kpc 
from the Sun, is one of the best laboratories to investigate the 
process of massive star formation. Several observational cam- 
paigns were performed with single-dish antennas and interfer- 
ometers (Codella et al. 1997; Furuya et al. 2002; Cesaroni et 
al. 2003; Belti-an 2004, 2005, 2006, 2007, 2011; MoscadelH et 
al. 2007; Vig et al. 2008) toward this region. G24 is associ- 
ated with a cluster of high-mass YSOs in different evolutionary 
stages, distributed in a region with size ~ 10" (see YSOs labelled 



Al, A2, B, C, and D in Fig. 1). Beltran et al. (2005, 2011) re- 
solved Al and A2 into five cores (see labels in Fig. 2) aligned 
in a southeast-northwest direction, which suggests a preferential 
direction for star formation in this region. The A cores are em- 
bedded in huge (0.1 pc), massive (a few 100 Mq) toroids, that 
rotate around the southeast-northwest (SE-NW) direction, and 
could still host elusive accretion disks in their interior (Beltran et 
al. 2005). This possibility would be even more intriguing for the 
object labelled Al, which is unique, because of the simultaneous 
presence of almost all the "ingredients" expected in a typical star 
formation "recipe": a 20 Mq star surrounded by a hypercompact 
Hii region, that is located at the centre of a rotating toroid un- 
dergoing infall towards the star The association of the rotation 
of toroids with the coherent motion of accretion disks would be 
supported by the detection of outflows, possibly driven by jets, 
along the SE-NW direction. The occurrence of outflows towards 
A and C has been first reported by Furuya et al. (2002) and has 
recently been imaged by Beltran et al. (2011) at high-angular 
resolution using CO emission. In particular, Beltran et al. (201 1) 
suggested that A2, and not Al, is the driving source of a bipo- 
lar CO outflow located along the SE-NW direction. However, 
given the complex outflow structure towards cores Al and A2, 
we cannot entirely discard the possibility that the core Al could 
be powering an additional outflow in the region. 

This doubt calls for high-angular resolution observations of 
a reliable jet tracer Silicon monoxide (SiO) thermal emission 
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is the best tool for this purpose: unhke other species such as 
CO, it is associated with shocks inside jets, suffers minimal con- 
tamination from infalling envelopes or swept-up cavities, and 
traces regions close to protostars, which are heavily extincted 
even in the near- and mid-infrared regions (e.g. Codella et al. 
2007; Lee et al. 2007a). Only SiO line emission will allow us to 
probe the mass loss process and unambiguously determine the 
direction of the flows driven by the G24 cluster. The formation 
of SiO is attributed to the sputtering of Si atoms from charged 
grains in a magnetised C-shock with velocities higher than 20 
km S-' (Schilke et al. 1997; Gusdorf et al. 2008ab). Although 
high-angular resolution studies of SiO in high-mass star-forming 
regions still refer to a quite limited number of objects (Hunter et 
al. 1999; Cesaroni et al. 1999; Qiu et al. 2007; Zhang et al. 2007), 
they confirm the power of SiO in tracing the mass loss process 
in complex environments like those typical of the massive star- 
forming regions. 

In this paper, we present SiO(l-O) and SiO(5^) im- 
ages obtained with the NRAOH Very Large Array (VLA) and 
SubMillimeter Array (SMAfl as well as a complementary 
30-m observations to unveil the mass loss process driven 
by the G24 cluster of high-mass YSOs. 

2. Observations 

2.1. VLA 

The G24 cluster was observed with 27 antennas of the NRAO 
VLA to measure the SiO(l-O) line emission at 43.4 GHz as well 
as the continuum emission. The observations were carried out in 
the Q-band with the D-configuration on August 10, 2008, and 
October 10, 2009, and with the C-configuration on August 14, 
2009. The half power beam width (HPBW) of the antennas is ~ 
1', which is the field-of-view of the images. The largest structure 
visible in the C+D-configuration is ~ 43". The phase reference 
centre of the observations was set to aaooo =18*^ 36" 12'5660, 
(52000 =-07° \2' W:\5. 

The SiO(l-O) line and the continuum were observed si- 
multaneously with the correlator in 2 IF mode. The continuum 
emission was collected with a 200 MHz bandwidth centred at 
43.339 GHz. The SiO(l-O) line was covered by simultaneously 
using two slightly overlapping 6.25 MHz (43 km s ') band^ 
Unfortunately, technical problems did not allow us to properly 
trace the SiO(l-O) emission in the VLA-D configuration. The 
observations were performed in fast-switching mode. Bandpass 
and phase were calibrated by observing 1 832-105, while the flux 
density scale was derived by observing 1331-1-305 (3C286). All 
data editing and calibration were carried out using the NRAO 
AIPSQpack age. The line cubes were obtained by subtracting the 
continuum from the line. Images were produced using natural 
weighting: Details of the synthesised cleaned beam, spectral res- 
olution, and rms noise of the maps are given in Table 1 . 

' The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc 

^ The SubMillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics, and is funded by the Smithsonian 
Institution and the Academia Sinica. 

3 IRAM is supported by INSU/CNRS (France), MPG (Germany), 
and IGN (Spain). 

* The 12.5 MHz bandwidth, needed to cover the whole SiO( 1-0) pro- 
file, was not available during the observations. 
' http ://w w w. aips.nrao.edu/index. shtml 



2.2. SMA 

The target was observed in the SiO(5^) line at 217. 1 GHz and in 
the continuum at 1 .4 mm with the SMA using two different array 
configurations. Compact- and extended-array observations were 
taken on May 21 and July 19, 2008, respectively. The covered 
baselines have lengths between ~ 8 and 100 k/i (compact) and 17 
and 160 k/1 (extended). We used two spectral sidebands, both 2 
GHz wide, separated by 10 GHz, covering the frequency ranges 
of 215.4-217.4 and 225.4-227.4 GHz, with a uniform spectral 
resolution of about 0.5 km s"' . The phase reference centre of the 
observations was set to 0-2000 = 18'' 36™ 12'.565, ^2000 = -07° 12' 
^.'90. 

Absolute flux calibration was derived from observations of 
Titan and Uranus. The bandpass of the receiver was calibrated 
by observations of the quasars 3C279 and 3C454.3. Amplitude 
and phase calibrations were achieved by monitoring 1743-038 
and 1911-201. We estimated the flux-scale uncertainty to be 
better than 15%. The visibilities were calibrated with the IDL 
superset Additional imaging and analysis was performed 

with MIRIAD (Sault et al. 1995) and GILDASQ . The continuum 
was constructed in the (u,v)-domain from the line-free channels 
both in the LSB and USB. Continuum maps were created by 
combining the data of both compact and extended configurations 
with the ROBUST parameter of Briggs (1995) set equal to zero, 
whereas SiO(5^) line channel maps were created using natural 
weighting after continuum subtraction in the UV plane. Angular 
and spectral resolution and as map sensitivity are given in Table 
1. 

2.3. IRAM 30-m 

Single-dish observations to prepare the interferometric observa- 
tional campaign were obtained with the IRAM 30-m telescope at 
Pico Veleta (Granada, Spain). The observations were carried out 
on January 27, 2007 pointing the telescope towards the position 
used as phase centre for the SMA observations: ffaooo = 18'' 36™ 
1Z565, (52000 = -07° 12' ia.'90. The pointing was checked by 
observing nearby planets or continuum sources and was found 
to be accurate to within 3-4". The observations were made by 
position-switching in wobbler mode. As spectrometer, an auto- 
correlator split into different parts was used to allow simultane- 
ous observations of four lines: SiO(2-l) at 86.8, (3-2) at 130.3, 
(5^) at 217.1, and (6-5) at 260.5 GHz, respectively (see Fig. 
3). The angular (HPBW) and the velocity resolutions provided 
by the backend, and the reached sensitivities are shown in Table 
1 . The integration time (ON+OFF source) was 70 minutes, while 
the main beam efficiency varies from about 0.77 (at 87 GHz) to 
0.48 (at 260 GHz). The spectra were calibrated with the standard 
chopper wheel method (uncertainty ~ 10%) and are reported 
here in units of main-beam temperature (T^mb)- 

3. Results 

3.1. Continuum emission 

Figure 1 shows the VLA and SMA maps of the continuum emis- 
sion at 7 and 1 .4 mm, respectively, towards G24. Table 2 sum- 
marises the position and the peak intensity of the detected cores. 
Given the lower angular-resolution, the present continuum im- 
age at 1 .4 mm does not add information with respect to our pre- 
vious PdBI and SMA observations at similar wavelengths (1.3- 

* http://cfa-www.harvard.edu/ cqi/mircook.html 
' http://www.iram.fr/IRAMFR/GILDAS 
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Table 1. List of the continuum and SiO observations 



Observation 


Telescope 


v" 






HPBW 


PA 


Spectral resolution 


1, 

rms noise 








(GHz) 


(K) 


(D-) 


(arcsec) 


(deg) 


(km s ') 


(mJy beam"') 


(mK) 


continuum 


VLA-C+D 


43.339 






1.5 X 1.1 


-11 




0.8 


309 


continuum 


SMA 


219.601 






1.5 X 1.4 


74 




6.0 


72 


SiO(l--O) 


VLA-D' 


43.424 


2 


9.6 


2.2 X 1.7'- 


-15^ 


0.67 


5 


870 


SiO(2-l) 


lRAM-30m 


86.847 


6 


19.3 


28 




0.54 


77 


16 


SiO(3~2) 


lRAM-30m 


130.269 


13 


28.9 


19 




0.36 


90 


18 


SiO(5^) 


lRAM-30m 


217.105 


31 


48.1 


11 




0.43 


173 


37 


SiO(5-4) 


SMA 


217.105 


31 


48.1 


1.7 X 1.4 


67 


0.67 


36 


415 


SiO(6-5) 


lRAM-30m 


260.518 


44 


57.7 


9 




0.36 


189 


42 



" Frequencies and spectroscopic parameters of the molecular transitions have been extracted from the Jet Propulsion Laboratory molecular 
database (Pickett et al. 19981. * For the molecular line observations the la noise is given per channel. The SiO(l-O) emission has been 
successfully observed only with the VLA-D configuration. 



Table 2. Position and peak intensity of the continuum cores" 



Core 


a(J2000) 


<5(J2000) 


pyeak ^peak 
7mm 1 .4mm 




(h m s) 


C ' ") 


(mJy beam" ) 


G24 Al 


18 36 12.55 


-07 12 11.05 


69 422 


G24B 


18 36 12.65 


-07 12 15.02 


12 


G24C 


18 36 13.11 


-07 12 07.42 


95 


G24D 


18 36 12.17 


-07 12 06.15 


22 



" The positions are based on the SMA image for all sources but B, which 
is detected only with the VLA. 

Table 3. Positions of the SiO emission peaks in the VLA and 
SMA images 



Peak 


a(J2000) 


<5(J2000) 




(h m s) 


(° ' ") 




SiO(l-O)- 


VLA 


Red 


18 36 12.55 


-07 12 10.54 


Blue" 


18 36 12.52 


-07 12 09.12 




SiO(5^) - 


SMA 


Red 


18 36 12.56 


-07 12 10.50 


Blue 


18 36 12.45 


-07 12 10.30 



" The SiO(l-O) blue-shifted emission is affected by absorption at low- 
velocities (see text). 

Table 4. List of the species and transitions serendipitously de- 
tected 



Transition 










(GHz) 


(K) 


(D^) 


C2H5CN(10i,,o-9,,5) 


86.820 


24 


0.5 


(CH3)2CO(184.14-175,13)-AE 


215.881 


111 


655.7 


'3CH30H-E(42 2-3l2) 


215.887 


45 


3.5 


CH3OCHO-E(20i, 20-19, i,) 


215.892 


298 


53.0 


CH30CHO-A(192.i8-18i.i7) 


216.360 


109 


7.2 


'3CH3OH-A(10_2.9-9-3.6) 


216.370 


162 


2.7 


H2 S (22,0-2 1, i) 


216.710 


84 


2.0 


C2H3CN(232,22-222,2,) 


216.937 


134 


996.4 


CH30H-E(5-i,4^-22) 


216.946 


56 


1.1 


CH30CHO-E(173,i4-163 n) 


216.959 


286 


43.7 


CH3OCHO-E(20i,20-19i 19) 


216.965 


111 


52.8 


CH3OCHO-A(20i.20-19i.i9) 


216.966 


111 


52.8 


'^^CH3OH-A(102.8-93.7) 


217.400 


162 


2.7 


HCC'''CN(24-23) 


217.420 


130 


346.4 


C2H5CN(295,25-285,24) 


260.536 


215 


417.0 



" Frequencies and spectroscopic parameters of the molecular transitions 
have been extracted from the Jet Propulsion Laboratory molecular 
database (Pickett et al. 19981 for all transitions except those of 
methanol, which have been extracted from the Cologne Database for 
Molecular Spectroscopy (Miiller et al. 12005b . 



1.4 mm; Beltran et al. 2005, 2011), which are consistent with 
the new observations. On the other hand, the present VLA beam 
is slightly better than that of our previous observations at 7 mm 
(2'.'3xr.'7, Furuya et al. 2002), but the sensitivity is a factor 2 
lower Again, the present VLA results agree with the previously 
performed analysis. The source B has been detected at 7 mm and 
not at 1 .4 mm is consistent with its evolved stage, which is asso- 
ciated with an ultra compact (UC) Hii region (see e.g. Furuya et 
al. 2002). But the analysis of the continuum emission is beyond 
the scope of the present paper, and will not be further pursued. 



3.2. SiO outflows 

The outflow activity was previously detected using emission of 
CO isotopologues by Furuya et al. (2002) and Beltran et al. 
(201 1) with the PdBI and SMA interferometers. Two molecular 
outflows, both oriented in the northwest-southest direction, were 
imaged. One (hereafter called outflow A) associated with the re- 
gion of the Al -I-A2 cores and is probably driven by A2. Another 
bipolar outflow is associated with core C. Figure 2 shows the 
maps of the integrated blue- and red-shifted SiO(l-O) emission 
observed with the VLA in the D configuration towards the G24 
cluster, while Fig. 4 shows examples of the SiO(l-O) profiles 
observed towards the positions of Al and of the blue- and red- 
shifted SiO(l-O) emission peaks of the outflow A (see Table 3). 
The cloud systemic velocity is -i-l 1 1 km s"', according to Furuya 
et al. (2002). The SiO(l-O) lines are characterised by emission 
up to high velocities (~ 20-25 km s"' with respect to the cloud 
velocity) in the wings, i.e. velocities comparable with what is 
observed using CO isotopologues. For outflow A, the VLA data 
confirm a poorly collimated and extended blue-shifted NW lobe 
and a smaller and more collimated red-shifted SE lobe. As seen 
in CO (Beltran et al. 2011), each lobe is characterised by the 
presence of a weaker extended counter-lobe. This could be a 
geometry effect because outflow A's main axis lies close to the 
plane of the sky. Alternatively, this could reflect the occurrence 
of a second outflow. For outflow C, the SiO(l-O) map shows a 
very elongated red-shifted (~ 25 km s"' w.rt. the systemic ve- 
locity) lobe: by dividing the length of the outflow by its width, 
we derive a coUimation factor fc ~ 4. In this case, we detect no 
SiO blue-shifted emission, whereas the CO data of Beltran et al. 
(201 1) clearly show blue lobe, albeit smaller (by a factor 3) than 
the red lobe. This suggests that outflow C could be located close 
to the edge of the molecular cloud and therefore the blue-shifted 
gas could flow through a low-density region. This could also ex- 
plain the lack of SiO emission, which is expected to trace the 
high-density collimated wind from the YSO. 
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G24.78+0.08, VLA, 7mm G24.78+0.08, SMA, 1.4mm 
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Fig. 1. Contour plots of the continuum emission at 7 (left panel) and 1.4 mm (right panel) imaged with the VLA and SMA inter- 
ferometers, respectively. The black triangles mark the positions of the sources identified by Furuya et al. (2002) and Beltran et al. 
(2004). The Icr rms of the maps is 6.0 (at 1.4 mm) and 0.8 mJy beam (at 7 mm), while first contours and steps correspond to 3cr 
and 6cr. The elhpses in the bottom-left corners show the HPBW: 1'.'5 x 1'.'4 (PA = 74°; SMA) and I'.'S x I'.'l (PA = -11°; VLA). 



Finally, no outflow activity has ever been observed towards 
source B, in agreement with its evolved stage, which is associ- 
ated with an UCHii region (Codella et al. 1997; Furuya et al. 
2002). No SiO emission (or CO emission) has been detected to- 
wards core D, which remains the most enigmatic object, bevause 
it is traced only by ~l-3 mm continuum emission and is not de- 
tected in any molecule (Furuya et al. 2002; Beltran et al. 201 1). 
Our non-detection seems to confirm that D is a non-centrally- 
peaked core without active star formation. 

The VLA image of the low-excitation (£„ = 2 K) SiO(l- 
0) emission is affected by absorption feature^ which prevents 
us from drawing a definite picture of the morphology of the 
SiO outflow and from identifying the driving source. Figure 5 
shows a zoom-in of the central region of G24. 78-1-0.08, cen- 
tred on the cluster of the A1-A2 continuum sources. The green 
contours correspond to the narrow absorption blue-shifted fea- 
ture (from -1-106 to -i-llO km s"', see Fig. 4) observed towards 
Al, against the hypercompact Hii region imaged by Beltran et 
al. (2007). This is ionised by a zero-age main-sequence star of 
spectral type 09.5 and mass of about 20 Mq. The hypercom- 
pact Hii region is very bright and therefore it is expected to ob- 
serve colder molecular gas in absorption against it. Indeed, at 
cm-wavelengths, where the brightness temperature of Al is ~ 
5 10^ K in a ~ a.'8 beam (Beltran et al. 2007), the NH3(2,2) 
spectral pattern has allowed us to trace the blue-shifted outflow 
lobe and the red-shifted gas infalling towards Al in absorption 
(Beltran et al. 2006). In conclusion, the VLA SiO(l-O) map cafls 
for the spectral analysis of a higher-J SiO at higher frequencies. 



* Unfortunately, we lack the information on the SiO velocity between 
+ 111 and +113 km s"' towards the position of Al because of technical 
problems. 



than those observed at SMA, where the brightness of the Al hy- 
percompact Hi: can be neglected. 

Figure 6 shows the SMA maps of the integrated blue- and 
red-shifted SiO(5^) emission. Only a hint of the outflow C is 
seen, and, as in the case of SiO(l-O), no emission has been 
observed towards cores B and D. On the other hand, the bipo- 
lar outflow A is clearly imaged, with the emission peaks lying 
closer to the A1+A2 system. The SiO(5^) structure is consis- 
tent with that traced by the SiO(l-O) line. Table 3 reports the 
positions of the SiO emission peaks detected in the VLA and 
SMA images. The position of the red-shifted peak is basically 
the same (within the uncertainties) for the for the SiO(l-O) and 
(5-4) lines, whereas, given the absorption observed at low veloc- 
ities that affects the SiO(l-O) map (see Fig. 4), the positions of 
the blue-shifted peaks in the two lines differ by ~ l'.'6. Example 
of SiO(5^) spectra are shown in Fig. 4. 

Figure 7 shows the maps of the SiO(5^)/SiO(l-0) inten- 
sity ratio, obtained after smoothing the maps to the same angu- 
lar resolution, and derived where both emissions have a signal- 
to-noise (S/N) > 3. The velocities where SiO(l-O) shows ab- 
sorption features (see Fig. 4) have not been taken into account. 
The red-shifted emission does not show a clear trend, while the 
SiO(5-4)/SiO(l-0) ratio for the blue-shifted emission peaks in 
correspondence of A2, suggesting an increase of the SiO excita- 
tion conditions. Finally, in Fig. 8 we compare the SiO(5^) pro- 
file observed with the IRAM 30-m antenna, with that obtained 
by integrating the SMA image over the IRAM beam (11"; blue 
histogram). The spectral shapes well match and, given the flux- 
scale uncertainties of the SMA (15%) and IRAM 30-m (10%), 
the line intensities also agree well, indicating that the SMA ar- 
ray recovers at least the 90% of the emission detected with the 
single-dish. 
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Fig. 2. Contour map of blue- and red-shifted SiO(l-O) VLA emission superimposed on the 1.3 mm continuum emission as observed 
at high angular resolution by Beltran et al. (2011) using a very extended SMA configuration ((X'SS x 0'.'44). The sources of the 
G24. 78-1-0. 08 cluster are labelled following Furuya et al. (2002) and Beltran et al. (2004, 2011). The SiO emission was averaged 
over the velocity intervals (h-80,+1 1 1) km s"' and (h-1 1 1,H-132) km s"' for the blue- and red-shifted emission, respectively. The Icr 
rms of the SiO maps is 1.2 mJy beam"' km s"'. Contour levels range from 5cr by steps of 3cr. The filled ellipse in the bottom-right 
corner shows the synthesised beam (HPBW): 2'.'2 x 1'.'7 (PA = -15°). 
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Fig. 4. Examples of SiO(l-O) and (5-4) line profiles (in brightness temperature Tb, scale) observed towards Al (left panel) and the 
peak positions of the blue- and red-shifted (middle and right panel) SiO( 1-0) line emission. Dashed lines mark the systemic velocity 
(H-111 kms-i). 



3.3. The origin of tine outfiow A 

Which is the driving source of the outflow A? The VLA and 
SMA SiO images have been overlaid on the 1.3 mm continuum 
map obtained at sub-arcsecond angular resolution (0'.'55 x 0'.'44) 
by Beltran et al. (201 1), which gives the best picture so far avail- 
able of the YSO population in the G24.78h-0.08 cluster It is 



possible to see that cores Al and A2 are resolved into three 
(called Al, Alb, Ale) and two (A2, A2b) cores, respectively. 
By analysing the CO emission, Beltran et al. (2011) suggested 
that the outflow A is actually powered by core A2, which is mas- 
sive (22 Mq) and associated with a hot-core (~ 180 K) detected 
in CH3CN, a typical hot-core tracer. The SiO images suggest 
that Al is not the driving source of the outflow motions detected 
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Fig. 6. Contour map of blue- and red-shifted SiO(5^) SMA emission superimposed on the 1.3 mm continuum emission as observed 
at high angular resolution by Beltran et al. (2011) using a very extended SMA configuration ((X'SS x 0'.'44). The Right panel reports 
a zoom-in of the central region. The sources of the G24. 78-1-0. 08 cluster are labeled following Furuya et al. (2002) and Beltran et 
al. (2004, 2011). The SiO emission was averaged over the velocity intervals (-(-89,-1-111) km s"' and (-1-111,-1-130) km s"' for the 
blue- and red-shifted emission, respectively. The rms Icr of the SiO maps is 18 mJy beam ' km s Contour levels range from 5cr 
by steps of 3cr. The filled ellipse in the bottom-right corner shows the synthesised beam (HPBW): 1'.'7 x I'.'A (PA = 67°). 



so far. Given the presence of several YSOs in the region, we 
cannot exclude additional fainter outflows. However, the peaks 
of the SiO(5^) emission suggest that the geometrical centre of 
outflow A, and thus its driving source, is located in the south- 
ern portion of the A2 core, which indeed is clearly elongated 
and probably hosts multiple YSOs. Additonal evidence comes 
from the position-velocity (PV) diagram of the SiO(5^) tran- 
sition along the NW-SE axis (see Fig. 9), which seems to con- 
firm that Al cannot be the driving source. Interestingly, the av- 
erage velocity of the SiO emission increases as a function of 
distance from the geometrical centre, suggesting gas accelera- 
tion. Figure 10 compares the PV diagram of SiO(5^) with that 
derived from the CO(2-l) SMA data by Beltran et al. (2011). 
We can clearly see that the CO emission (i) is completely ab- 
sorbed by foreground material at low velocities and is affected 
by filtered extended emission, confirming the need of a tracer 
such as SiO to properly image the molecular outflow, and (ii) 
shows additional components with respect to what was traced 
by SiO. Indeed, bright CO emission is seen at (-5", -1-115 km 
s~'), not traced by SiO, as well as weak emission at (-i-2",-h113 
km s ' ). In other words, the CO emission cannot be simply inter- 
preted as an intensity-scaled version of the SiO emission. These 
findings are not surprising given what already found for the jets 
driven by Sun-like protostars, where, thanks to the smaller spa- 
tial scales that can be investigated, it is possible to see that CO 
traces not only the jet but also the walls of the cavity opened by 
the jet itself (e.g. Lee et al. 2007b, see their Fig. 5). In summary, 
the present SiO maps support a formation mechanism for mas- 
sive stars similar to that of their low-mass counterparts, where 
jets, clearly traced by SiO emission, create outflows of swept- 
up ambient gas traced by CO. The possibility of having an SiO 
jet is also supported by the extremely high velocity (up to ~ 70 
km s"' with respect to the systemic velocity) of the SiO emis- 
sion detected in the 7 = 2-1 and 3-2 spectra thanks to the high 



sensitivities provided at these frequencies by the 30-m antenna. 
Of course, in the present case, the angular resolution is not high 
enough to assess the occurrence of SiO jets. However, in anal- 
ogy to the SiO jets in low-mass stars, one cannot rule out the 
possibility that the elongated SiO outflow does trace jet activity. 



3.4. Other molecular species: the outflow and the hot core 

The 2-GHz-wide LSB bandwidth used to trace SiO(5^) emis- 
sion with the SMA allowed us to serendipitously observe several 
lines of different molecular species, which are listed in Table 4 
and shown in Fig. 11. In particular, the LSB spectrum is dom- 
inated by broad H2S(22.()-2i.i) emission, which confirms that 
hydrogen sulphide is an excellent tracer of molecular outflows 
(Codella et al. 2003, Gibb et al. 2004). Figure 12 shows the inte- 
grated blue- and red-shifted H2S emission from outflow A, while 
Fig. 9 shows the corresponding PV diagrams. 

Outflow A is clearly seen in the H2S image, which sug- 
gests (like SiO) that the geometrical centre (and thus the driv- 
ing source) is located in the southern portion of the A2 core. 
Interestingly, the H2S outflow is definitely less collimated than 
the SiO(5^) one, observed with the same angular resolution. 
H2S and SiO show different behaviours also in the PV diagrams, 
where hydrogen sulphide does not show a clear increase of the 
average velocity with the distance from the driving source. These 
findings could be the signature of an additional fainter outflow 
driven by one of the A1-I-A2 YSOs that contributes to the ob- 
served emission. However, what we found agrees with an en- 
hancement of the H2S abundance as a consequence of the evap- 
oration of the dust mantles in a shocked gas (e.g. van Dishoeck 
& Blake 1998), whereas SiO comes from a definitely smaller re- 
gion that is directly associated with the primary jet, where the 
refractory dust cores are disrupted as well. 
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Fig. 7. SiO(5^)/SiO(l-0) intensity ratio (brightness temperature, T^, scale), derived where both emissions have an S/N > 3. For this 
comparison, the SMA image has been smoothed to the angular resolution of the VLA image. The velocities affected by SiO(l-O) 
absorption (see Fig. 4) have not been taken into account: right and left panels are for the blue (+89,+ 105 km s"') and red (+1 15,+ 130 
km s"') emission, respectively. Contour levels range from 0.5 to 2.0 by steps of 0.5. The sources of the G24.78+0.08 cluster are 
labelled following Furuya et al. (2002) and Beltran et al. (2004, 201 1). 



In addition to SiO and H2S, we detected several high- 
excitation {Ea between 45 and 162 K) transitions of methanol 
isotopologues (see Table 4). Figure 11 shows the spectra ob- 
served towards the A2 position. The CH3OH line profiles are 
narrower than those of SiO and H2S and, although a weak com- 
ponent related to outflows cannot be excluded for '^CHsOH, it 
appears that methanol emission is dominated here by gas heated 
by YSOs. Indeed, from Fig. 9 one clearly sees that the PV di- 
agram of methanol, outlining a compact circular pattern, is sig- 
nificandy different from that of SiO. Most likely, the CH3OH 
molecules trace the hot core A2 (already imaged by Beltran et 
al.(2005), and are released from dust mantles because of stel- 
lar irradiation. In particular, '^CHsOH peaks at a position that 
matches the geometrical centre of the outflow well (see the solid 
vertical line in Fig. 9). Thus, the present observations suggest 
that a YSO lying in the SB of core A2 is heating the gas and 
driving the outflow A. 



Finally, we also detected several lines in the SMA spectral 
window that are due to cyanoacetylene (in HCC'^CN form) as 
well as to other complex molecular species usually considered as 
hot-core ti-acers: CH3OCHO (methyl formate), C2H3CN (vinyl 
cyanide), and, tentatively, one transition of (CH3)2CO (acetone). 
Other complex species were reported by Beltran et al. (2011, see 
their Table 1). Indeed, these emissions are characterised by high 
excitation (E^ in the 109-298 K range) and compact spatial dis- 
tribution (see the C2H3CN in Fig. 9 as an example), indicating an 
association with the hot core traced by '^CH30H. Remarkable 
is the tentative detection of acetone, given that only recently 
this organic specie has been detected towards the hot molecu- 
lar core Sagittarius B2(N-LMH) and in Orion BN/KL (Snyder et 
al. 2002, Friedel et al. 2005, Goddi et al. 2009). 



4. Physical conditions as traced by SiO 

4.1. SiO(1-0) absorption 

From the SiO(l-O) absorption feature one can obtain an esti- 
mate of the excitation temperature of this transition. We assume 
that the absorption is due to the blue lobe of the outflow lying 
between the observer and the hypercompact Hii region. Since 
the lobe is resolved in our maps, while the Hii region is much 
smaller than the synthesised beam of the SiO images, the SiO(l- 
0) brightness temperature measured in the synthesised beam is 
a mixture of absorption (towards the HII region) and emission 
(from the rest of the beam): 



nerB = (nB-ns)re + Qsr„ 



(1) 



where and are the brightness temperatures of the SiO(l-O) 
line in emission and absorption, while Qs and Qb are the solid 
angle subtended by the Hii region and the synthesised beam of 
the SiO(l-O) map. Moreover, Tg and are given by 



T, = (rex-rHii)(i -e-^^). 



(2) 
(3) 



with Tex excitation temperature of the 7=1-0 transition and Tun 
brightness temperature of the continuum emission from the hy- 
percompact Hii region. Here, we have assumed that Tex is con- 
stant along the line of sight and across the region covered by the 
synthesised beam. After some algebra, one obtains the expres- 
sion 



1 



Qb 1 - e 



(4) 
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Fig. 9. Position-velocity cut of SiO(5^), H2S(22,o-2i,i), CH30H-E(5-i.4-4_22), and i3CH3OH-A(102,8-93,7) along the outflow A 
(PA = -40°). The position offsets are measured from the Al position, positive towards northwest. The rms Icr is 20 mJy beam"'. 
Contour levels range from 3cr by steps of 3cr for SiO and '^CHsOH-A, and from 3cr by steps of lOcr for H2S and CH3OH-E. 
Dashed vertical lines mark the position of the Al and A2 YSOs, while the solid vertical line is for the average peak position of the 
typical hot-core tracers (see text). Solid horizontal line shows the systemic velocity (-1-III.O km s '). 



We note that the Qb is relatively small compared to the size of 
the outflow lobe, and hence the opacity should not depend signif- 
icantly on the line of sight. Also, for obvious reasons, the Hii re- 
gion either lies behind the lobe or is enshrouded by it. Therefore, 

Qs 
Qb 

written as 



< 1. Since ^p- <s; 1, Eq. dUl can be 



' HII 



(1-e-a). 



(5) 



Most likely, the absorption and emission lines of sight are cross- 
ing similar amounts of gas, and we hence make the additional 
assumption that ^ t^. By replacing Eq. (|2]i into Eq. (|5]l, one 
can finally obtain an expression for T^^: 



(6) 



From our observations we measure absorption with Tb - 
-6 K in a beam of Q.^ - 4.24 arcsec^, while from the data 
of Beltran et al. (2007) one obtains Qs - 0.0554 arcsec^ and 
a source-averaged brightness temperature of the Hii region at 
7 mm Thii - 960 K. An estimate of T^. can be obtained from the 



SiO(l-O) spectra one beam away from the absorption peak: Tg - 
3-6 K. With these numbers we obtain Tgx - 4-8 K. The bright- 
ness temperature of the continuum emission measured at 7 mm 
towards the Hn region in our images is ~12 K, i.e. comparable 

to the value ( j^Tnu ^ 12 K) estimated applying beam dilu- 
tion to the high-resolution measurement of Beltran et al. (2007): 
this proves that basically the entire continuum emission in our 
beam is due to free-free emission from the ionised gas, while 
dust emission is negligible. 



4.2. LVG analysis 

We ran the RADEX non-LTE model (van der Tak et al. 2007) 
with the rate coeflicients for collisions with H2 reported by 
Turner et al. (1992) using the escape probability method for a 
plane-parallel geometry to fit the observed SiO line ratios and 
brightness. We explored H2 densities from 10^ to 10^ cm"^, ki- 
netic temperatures Tkin from 50 to 500 K, and an LVG optical 
depth parameter n(SiO)/ (dV/dz) - A^sio/AV ranging from 10'^ 
to 10'^ cm"- (km s"')"', i.e. from the fully optically thin to the 
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Fig. 10. Comparison between the position-velocity cut of SiO(5^), black contours, and that of CO(2-l) in colours, derived along 
the outflow A (PA = ^0°) using the CO dataset presented by Beltran et al. (2011). The solid horizontal line shows the systemic 
velocity (+111.0 km s"'). 



optically thick regime. We used an FWHM linewidth of 10 km 
s"', as suggested by the SiO spectra (see Fig. 4). 

As a first step, we modelled the SiO emission observed with 
the IRAM 30-m telescope, using the SiO(2-l) and (5^) tran- 
sitions, which are sensitive to similar excitation conditions (see 
Table 1) as the lines observed with the interferometers (discussed 
below). The low-J SiO lines, as already learned from previous 
studies of low-mass protostellar systems (e.g. Cabrit et al. 2007 
and references therein), are not very sensitive to kinetic tem- 
perature. Figure 13 thus reports the solutions for the observed 
SiO(5^)/SiO(2-l) intensity ratio at the typical red- (+116 km 
s"') and blue-shifted (+106 km s"') emission in the nu^-N{SiO) 
plane and for two extreme temperatures: 50 and 500 K. The ob- 
served intensity ratio was corrected for beam dilution assuming 
an emitting source of 7", which is representative of the sizes 
of the outflows imaged with the VLA and SMA. To show how 
much the results depend on this assumption on emitting size, 
we also report the solutions for a smaller size (i.e. 3"), which is 
the typical size of the brightest SiO clumps in G24. The dashed 
contours take into account the uncertainties associated with the 
intensity ratio. The 30-m spectra do not provide any constraint 
on the kinetic temperature and column density, while densities 
of below 10^ cm"-^ can be inferred. 

On the other hand, we can obtain tight constraints on vol- 
ume and column densities by analysing the SiO emission ob- 
served with the VLA and SMA arrays, after degrading the SMA 
map to the angular resolution of the VLA image. Figures 14 and 
15 report the solutions for the observed SiO(l-O) and SiO(5- 
4) brightness temperatures at the same red- (+116 km s"') and 
blue-shifted (+106 km s"') velocities investigated using the 30- 
m spectra. We modelled the emission observed towards three 



positions: Al (Fig. 14), as well as where the J=l-0 emission 
is only weakly afifected by absorption, i.e. at the position of the 
red- and blue-shifted SiO(l-O) emission peaks (Fig. 15). As re- 
ported in Table 3, the positions of the red-shifted peak as traced 
by SiO(l-O) and (5^), hence once considered the angular res- 
olutions, agree well. We report the solutions found for the same 
kinetic temperatures used above to model the 30-m emission (50 
and 500 K). We plot in black in Fig. 14 the contours correspond- 
ing to the the excitation temperature of the SiO(l-O) line es- 
timated in Sect. 4.1 from the SiO(l-O) absorption feature ob- 
served towards Al. 

In practice, the present LVG plots indicate for the three posi- 
tions (i) volume densities between lO-' and 10^ cm""*, and (ii) 
well-constrained SiO column density, in the 0.5-1 10'^ cm"^ 
range. These solutions are consistent with what was found using 
the 30-m IRAM SiO(5-4) and (2-1) spectra. Nisini et al. (2007) 
performed an LVG analysis of SiO emission towards the nearby 
prototypical Class low-mass objects LI 148 and LI 157, using 
single-dish (IRAM, JCMT) data on angular scales of 10". In that 
case as well, no tight constraints have been obtained for the ki- 
netic temperature, whereas the volume densities were found to 
be between 10^ and 10^ cm"^. Moreover, Gusdorf et al. (2008a) 
modelled the SiO emission from L1157 using a C-shock code, 
which led to pre-shocked densities of 10''-10^ cm""*. Thus, al- 
though completely different spatial scales are involved, the vol- 
ume densities here inferred for G24 (10^-10^ cm"^) appear to be 
consistent with those derived for L1448 and LI 157. 

Assuming that the physical conditions reported above are 
representative of the whole SiO outflow (including where 
SiO(l-O) emission is affected by absorption), using a volume 
density of 10"^ cm"^, the inferred size of the lobes (7"), and the 
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Fig. 3. SiO line profiles (in main-beam temperature, Tmb, scale) 
observed with the IRAM 30-m antenna towards G24.78+0.08. 
The dashed lines stand for the systemic velocity (+111.0 km 
s"'). Two ethyl cyanide (C2H5CN) emission lines have been 
serendipitously detected at 86819.85 MHz (10i,io-9i,9; £^-24 
K; bottom panel), and at 260535.69 MHz (295,25-285,2^; £'u=215 
K; upper panel, blended with the SiO(6-5) blue wing). 



G24 distance (7.7 kpc), the corresponding outflow mass is ~ 40 
Mq. This implies, according to Beuther et al. (2002), Zhang et 
al. (2001, 2005), and Lopez-Sepulcre et al. (2009, see their Fig. 
5), a luminosity for the driving source of the SiO outflow of a 
few 10"^ Lq, coiTesponding to a late O-type ZAMS YSO (e.g. 
Panagia 1973). If we conservatively assume riu, = lO-' cm"-', 
the driving source is still a massive YSO, corresponding to a 
B2 spectral type. However, for consistency purposes, we cross- 
checked the outflow masses using the estimates of SiO column 
density and assuming a given SiO abundance relative to H2. The 
latter has a large uncertainty since SiO can be greatly enhanced 
in shocks. Assuming typical abundances of 10^^-10"^ (i.e. 4- 
5 orders of magnitude greater than that in dark clouds, 10 
Ziurys et al. 1989), we inferred volume densities between 10"^ 
and 10^ cm"-', and thus again masses (~ 40^00 Mq) in agree- 
ment with an outflow driven by a late O-type ZAMS YSO. 



5. Conclusions 

We conducted a multiline SiO survey towards the G24.78+0.08 
region, which is an excellent laboratory to study the process of 
high-mass star formation, because it is associated with YSOs 
in different evolutionary stages. After preliminary IRAM 30-m 
single-dish runs, we obtained high angular resolution images us- 
ing the VLA and SMA interferometers. The main results can be 
summarised as follows: 
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Fig. 5. Zoom-in of the central region of G24.78+0.08, where 
the cluster of A1-A2 continuum sources (grey scale) has been 
mapped at high angular resolution by Beltran et al. (201 1) us- 
ing a very extended SMA configuration (0'.'55 x 0'.'44). The spa- 
tial distribution of the blue-shifted (from +106 to +1 10 km s"') 
absorption (see Fig. 4) observed towards Al is reported by the 
green contours. Contour levels range from -lOcr to -3cr by steps 
of Icr (0.3 mJy beam"' km s"'). The filled ellipse in the bottom- 
right corner shows the synthesised beam (HPBW): 2'.'2 x 1'.'? (PA 
= -15°). 
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Fig. 8. Comparison between the SiO(5^) line profile (in flux 
density, Sy, scale) observed with the IRAM 30-m antenna 
(black) and the spectrum obtained from the SMA data (blue) by 
integrating over a region equal to the IRAM 30-m beam (11")- 
For this comparison the IRAM 30-m spectrum has been resam- 
pled to the spectral resolution of the SMA data (0.67 km s"'). 
The dashed lines stand for the systemic velocity (+111.0 km 



1. High-velocity SiO emission (up to 25 km s"' w.rt. the sys- 
temic velocity, +111 km s"') reveals two collimated outflows 
driven by the A2 and C millimeter continuum massive cores. 
On the other hand, no SiO emission has been detected to- 
wards more evolved young stellar objects associated with an 
UCHii region (core B) or driven by the hypercompact (core 
Al) Hii regions. Moreover core D shows no SiO emission. 
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Fig. 13. Analysis of the SiO red-shifted line emission in the G24.78+0.08 outflow A observed with the IRAM 30-m antenna. The 
solutions for the observed SiO(5^)/SiO(2-l) intensity ratio are shown in the nH2-^(SiO) plot for non-LTE (RADEX) plane-parallel 
models at the labeled kinetic temperatures. Solid contours are for the measured ratios at the typical red- (-1-116 km s~') and blue- 
shifted (-1-106 km s ') velocities, after correction for a beam dilution derived assuming a source size of 3" and 7" (see text). The 
SiO(5^)/SiO(2-l) ratios are 0.21 (3"; +116kms-^), 0.19(3"; H-106kms-i), 0.27 (7"; -h1 16 km s"'), and 0.30 (7"; +106 km s'^). 
Dashed contours take into account the uncertainties. 



confirming its quiescent nature, without any signature of star 
formation. 

2. The LVG analysis of the SiO emission reveals high-density 
gas (10^-10^ cm~^), with clearly constrained SiO column 
densities (~ 10'^ cm"^). The average velocity of the SiO 
emission increases as a function of distance from the driv- 
ing source, which suggests gas acceleration. Although the 
angular resolution is not high enough to demonstrate the oc- 
currence of SiO jets, if we assume the standard approach that 
SiO is tracing shocks inside jets, it is reasonable to associate 
the observed colUmated SiO structures with jet activity. 

3. The driving source of the A2 outflow (i) has an estimated lu- 
minosity of > 10^ Lq, which is typical of a late O-type star, 
and (ii) is located at the centre of a hot molecular core (traced 
in the present data set by emission from methyl formate, 
vinyl cyanide, cyanoacetilene, and acetone) that rotates on 
a plane perpendicular to the outflow main axis. 

4. To our knowledge, we obtained one of the first interfero- 
metric images of an SiO jet-like outflow from young > 10"^ 
Lq stars. High spatial resolution maps of SiO high-velocity 
emission driven from young > lO'' L© stars have have been 
so far obtained towards IRAS18264-1152, 18566-1-0408, 
20126+4104, and 23151+5912 (Cesaroni et al. 1999; Qiu et 
al. 2007; Zhang et al. 2007). 1RAS20126+4104 is probably 
the best example of a jet from a massive YSO so far ob- 
served; it is traced by SiO, H2, and [Fell] emission (Caratti 
o Garatti et al. 2008). To conclude, the present SiO observa- 
tions support the theory that O-type stars form according to 
a core accretion model, i.e. via a scaled-up picture typical of 
Sun-Uke star formation, where jets, well-traced by SiO emis- 
sion, create outflows of accumulated and accelerated ambi- 
ent gas that in turn is well traced by CO. 
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Fig. 14. Analysis of the SiO line blue- and red-shifted emissions towards Al as observed with VLA and SMA (convolved to 
the angular resolution of the VLA image). The solutions for the observed SiO(l-O) and SiO(5-4) brightness temperatures (solid 
lines) at the typical red- (-1-116 km s"') and blue-shifted (-1-106 km s"') velocities are shown in the nu,-N{SiO) plane for non- 
LTE (RADEX) plane-parallel models at the labelled kinetic temperatures. Dashed lines are for the uncertainties. The SiO(5^) 
brightness temperatures are 4.4 K (-1-116 km s"') and 4.3 K (+106 km s"'), while the SiO(l-O) brightness temperature is 4.0 K 
(h-1 16 km s '). Given the blue-shifted absorption observed towards Al, we plot in black the solutions corresponding to the derived 
SiO(l-O) excitation temperature (^ 6 K, see Sect. 4.1). 




Fig. 15. Analysis of the SiO line blue- and red-shifted emissions towards the corresponding emission peaks of the SiO(l-O) map 
(see text), as observed with VLA and SMA (convolved to the angular resolution of the VLA image). The solutions for the observed 
SiO(l-O) and SiO(5^) brightness temperatures at the typical red- (-1-116 km s"') and blue-shifted (-1-106 km s"') velocities are 
shown in the nH2-^(SiO) plane for non-LTE (RADEX) plane-parallel models at the labelled kinetic temperatures. The SiO(l-O) 
brightness temperatures are 4.7 K (-1-116 km s"') and 6.1 K (-1-106 km s"'), while the SiO(5^) brightness temperatures are 1.1 K 
(-1-116 km s ') and 2.9 K (-1-106 km s '). Dashed lines are for the uncertainties. 
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Fig. 11. SiO, H2S, CH3OH, and '^CH^OH line profiles (in 
brightness temperature, Ts, scale) observed towards the A2 hot 
core with the SMA. The dashed lines denote the systemic ve- 
locity (+111 km s"'). Several high-excitation lines due to typi- 
cal tracers of hot-core chemistry (methyl formate, acetone, vinyl 
cyanide, cyanoacetilene) have also been detected (see Table 4). 
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Fig. 12. Contour map of blue- and red-shifted H2S(22,o-2i,i) 
SMA emission superimposed on the 1.3 mm continuum emis- 
sion as observed at high angular resolution by Beltran et al. 
(201 1) using a very extended SMA configuration ((X.'SS x 0'.'44). 
The sources of the G24.78+0.08 cluster are labelled following 
Furuya et al. (2002) and Beltran et al. (2004, 2011). The H2S 
emission was averaged over the velocity intervals (-1-94,-1-111) 
km s"' and (+111, +124) km s"' for the blue- and red-shifted 
emission, respectively. The rms Icr is 18 mJy beam"' km s"'. 
Contour levels range from 5cr by steps of lOcr. The filled ellipse 
in the bottom-right corner shows the synthesised beam (HPBW): 
1'.'7 X 1'.'4 (PA = 67°). 
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